Tau energy losses at ultra-high energy: continuous versus stochastic treatment. 
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We study the energy losses of the tau lepton in matter through electromagnetic processes at 
ultra- high energy (UHE). We use both a stochastic and a continuous framework to treat these 
interactions and compare the flux of tau leptons propagated after some amount of matter. We 
discuss the accuracy of the approximation of continuous energy losses by studying the propagation 
in standard rock of taus with both mono-energetic and power law injection spectra. 



I. INTRODUCTION 

The observation of high-energy cosmic neutrinos pro- 
duced in distant astrophysical sites (or possibly by other, 
more exotic, mechanisms) has become one of the major 
challenges of astroparticle physics. In both astrophysi- 
cal and exotic models, substantial fluxes of electron and 
muon neutrinos are expected from the desintegration of 
charged pions (and kaons) produced in the interaction of 
accelerated particles with ambient matter and radiation, 
either at the source location or on their way through the 
universe. Given the large distances traveled by these cos- 
mic neutrinos, approximately equal fluxes in f,,,^ and 
v T are expected on Earth as a result of flavour mixing 
and oscillations [H, @, 0] • 

At ultra-high energies, the Earth becomes opaque to 
neutrinos as a result of charged- and neutral-current in- 
teractions which deplete their energy or convert them 
into the associated charged lepton in case of charged- 
current. Tau neutrinos can then be traced through the 
observation of byproducts of the decay of the associated 
tau lepton in matter such as ice 0, , water @, [1] , or 
in the air 0, [l(| • 

The tau lepton is subject to energy-dependent radia- 
tive processes : Bremsstrahlung, production of e + e~ 
pairs and photonuclear interactions (ionisation is neg- 
ligible at ultra-high energy). A comprehensive knowl- 
edge of the probability that a tau with initial energy Eq 
will not be absorbed during its propagation in matter is 
thus required to predict the expected rates of tau neu- 
trinos detection by various experiments. In the muon 
case, because of its large lifetime, the survival probabil- 
ity is entirely determined by the energy losses over the 
whole energy range. Neglecting the fluctuations in en- 
ergy losses (or, in other words, adopting the continuous 
energy loss approach) is known to lead to an overestimate 
of the survival probability In contrast, due to the 

short lifetime of the tau, the calculation of the range in 
that case is determined by the decay up to « 10 8 GeV. 
Then, the attenuation length becomes smaller than the 
decay length. Fluctuations in the energy losses may then 
play an important role - as in the muon case. 

In this respect, it has already been put forward that 
such effects may indeed not be negligible for the tau [12] • 



The approach adopted in this previous work was to in- 
ject tau neutrinos in matter as a source of tau leptons, 
and then to compute the flux of tau leptons after some 
fixed distance of matter. This calculation was performed 
in two ways: the first by using a full simulation of all 
involved processes through a Monte-Carlo generator to 
account for the stochasticity of the problem, and the sec- 
ond by using a simplified semi-analytical framework to 
reproduce the solution within the continuous energy loss 
framework. Eventual differences in the flux calculations 
were attributed to the lack of accuracy of the continuous 
energy loss framework. 

Such a strategy however makes it difficult to disentan- 
gle the direct study of t propagation from other effects 
related with the propagation and interactions of the v T . 
The aim of this paper is thus to re-examine the effects 
of the fluctuations in energy losses due to radiative pro- 
cesses of the tau lepton at ultra-high energy in a frame- 
work where we choose to directly inject r's instead of 
2/ T 's. This method allows to simplify the problem and 
to get a more direct comparison between the continuous 
and the stochastic treatment of the energy losses. 

In section UU we present the physical processes that we 
take into account in this study, and the methods we will 
use to compute the propagated flux within both the con- 
tinuous energy loss framework and the stochastic one. 
In section IIII1 we apply our calculations to the case of 
mono-energetic beams of tau leptons, whereas in section 
IIVI we do the same in case of a power law injection spec- 
tra of tau and compare the results obtained within the 
two frameworks. Finally, we discuss our conclusions in 
section [V] 



II. PROPAGATION OF TAU THROUGH 
MATTER 

A. The processes 

The tau lepton is an unstable particle. In the observer 
frame, its decay length is energy-dependent : 

\ dec (E)=49(-^\m. (1) 
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At the energies of interest for this study, the relevant 
electromagnetic processes that the tau lepton undergoes 
are the Bremsstrahlung, the pair production and the pho- 
tonuclear interactions. Whereas the cross sections of the 
two first ones are well established [TH [l4[ , the description 
of the photonuclear process relies on the proper modeli- 
sation of the nucleon structure functions at low x and 
high Q 2 , and several calculations of the corresponding 
cross section can be found in the literature. The dif- 
ferences reside in the parameterisation of the relevant 
structure functions obtained from different formalisms. 
A popular treatment uses a combination of the gener- 
alized vector dominance model for the non-perturbative 
regime [HI , and of the color dipole model for the pertur- 
bative part [IB]. Other widely used results are obtained 
on basis of parameterisations of data using Regge the- 
ory [13, (3 • A recent comparative study shows that all 
these calculations give rates of photonuclear energy loss 
in good agreement except for [2(| which results in 
a significantly higher rate at energies above 10 8 GeV. A 
new calculation based on saturation physics [2l| is also 
presented in [l9j . In this case the corresponding energy 
loss rate is more than a factor of 2 lower than the stan- 
dard results at the highest (~ 10 12 GeV) energies. In the 
following we use the parameterisation obtained by Dutta 
et al. [l7| for all the calculations we are presenting. 

Finally, the tau lepton is also subject to weak interac- 
tions. However, using the cross sections derived in [22T ]. 
the corresponding interaction length is expected to be 
much larger than for the electromagnetic processes un- 
til an energy greater than 10 12 GeV. As we are only 
interested in studying the effect of the fluctuations of 
the electromagnetic interactions of the tau, we can thus 
comfortably neglect its weak iteractions. We do as well 
regarding the regeneration of t leptons at lower energy 
through the chain r — > v r —> t via a tau decay followed 
by a charged-current weak interaction of the neutrino. 



B. Solving the equation of transport 
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FIG. 1: Evolution of f3 as a function of the energy for a 
t lepton in standard rock (A=22). From bottom to top: 
bremsstrahlung, pair production and photonuclear interac- 
tions. 



1. Continuous energy losses approximation 

When the differential cross-sections exhibit a peak near 
y=Q, as it is indeed the case for the processes we are deal- 
ing with, the integrals are dominated by the behavior of 
the integrands around 0, in such a way that an expansion 
of these integrands can be performed [23j . At first order 
in y, this yields the following equation : 

d<P(E,x) $(E,x) d ( „, , , \ 

where we have introduced the standard notation 



Within the above assumptions, the transport equation 
that describes the evolution of the t flux $ along its path 
trough matter, accounting for all processes of creation or 
absorption of a tau with an energy E at a position x, is 
given by: 
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where Adec(E) is the decay length, Ai(E) and 
d(Ji/dy{E 1 y) respectively the interaction length and the 
differential cross section of each process, y the inelastic- 
ity of the process considered, p the density of matter, M 
the Avogadro number and A the atomic mass number. 
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We show on Figfl] the evolution of j3 as a function of 
energy for the three processes described in the previous 
sub-section. Note that the contribution of the energy loss 
due to the photonuclear interactions does not reach any 
asymptotic value because of the expected growth of the 
photo- nucleon cross section <t 7 at. 

Within the approximation of continuous energy losses, 
the average energy lost per unit distance dE/dx is, at 
ultra-high energy, assumed to be proportional to the 
mean inelasticity of each process in the following way 



f x =-p^Ej2 { y i )a i (E)- 



(5) 
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FIG. 2: Comparison of the semi-analytical solution (solid line) 
of the transport equation within the continuous energy losses 
approximation (Eqn. [6} with respect to a Monte-Carlo calcu- 
lation (squared points) within the same approximation. An 
E~ 2 flux between 10 8 GeV and 3 ■ 10 11 GeV is used at in- 
jection. Three distances of propagation in rock are shown: 1 
km, 5 km and 10 km (from top to bottom). 



The right-hand side of this expression is nothing else but 
—pE(3(E). This implies that Eqn. [3] can be easily inte- 
grated, leading to the following expression : 

$(E,x) = <S>o(E ) exp [ X d u (JL 7 (E u ) - ) )(6) 

JO \9E \dec{E u )J 

where "/(E) = pE(3(E) and E v the solution of 
Je v l(E T ) 

In the following, we use Eqn. [5] to compute any propa- 
gated flux of tau leptons when referring to as the contin- 
uous energy losses approximation. To verify this expres- 
sion, we performed a Monte-Carlo calculation evaluating 
at each step the decay probability as well as the continu- 
ous energy losses through Eqn. [SJ For an incident flux of 
tau leptons following a E~ 2 spectrum in energy between 
10 8 GeV and 3 • 10 11 GeV, we show the agreement of the 
two calculations on Fig|5] after a propagation in 1, 5, and 
10 km of rock. Let us point out that, within the con- 
tinuous energy losses approximation, the sharp cutoffs 
in energy result from the choice of the maximal energy 
at x = 0, which is univocally related to the propagated 
energy at depth x. 



FIG. 3: Interaction lengths computed in the y range y € 
[ycut,y ma x] with y cut = 1CT 3 . 



2. Stochastic treatment 

To test the accuracy of results obtained by using the 
continuous energy losses approximation, we need to ac- 
count for the stochastic nature of the transport equation. 
A convenient way to solve Eqn. [5] is to use a Monte-Carlo 
generator sampling all the interactions. However, to limit 
CPU time spent where the cross sections are large but 
the energy losses are small, a standard way to proceed 
is to separate the losses into two components (TH, l24j : 
a continuous one where the rate of the losses is large 
(y £ [ymin, Ucut]), and a stochastic one where the differ- 
ential cross sections lead to more catastrophic losses but 
with a weaker rate (y 6 [y C ut, Umax})- In practice, this 
means that within an elementary step, the Monte-Carlo 
samples all the interactions according to cross sections 
computed for y > y cut together with the probability to 
decay. At the same time, a continuous energy loss is ap- 
plied according to Eqn. [U except that the upper bound 
of the (3 coefficient is replaced by the cut y cu t- In that 
way, the y range where the stochastic nature of the in- 
teractions can be relevant is taken into account, and at 
the same time, the frequent small losses due to the peak- 
ing of the cross sections near y — are applied. A good 
compromise to reproduce the stochastic features using a 
reasonably fast code is to take y cut = 10~ 3 . 

To get a feeling of the effects which may be induced 
by fluctuations in energy losses, we plot on Fig[3]the in- 
teraction length of each process, but restricting the y 
range to y G [y C ut,ymax]- The interaction length of the 
Bremsstrahlung is clearly the higher, so that the effects 
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induced by this process are marginal. On the other hand, 
the interaction lengths induced by pair production and 
photonuclear interactions are lower and rather close to 
each other, leading to higher rates. Stochastic effects are 
not expected to be strong for pair production processes, 
as the corresponding differential cross-section behaves as 
rs y~ 2 ■ On the contrary, the rate of photonuclear in- 
teractions is never negligible, and in particular, above « 
6 EcV, this process becomes the dominant one. Its dif- 
ferential cross section behaves as y~ 1+a with a slighty 
decreasing with energy from 0.1 at 10 s GeV to 0.05 at 
10 12 GeV. Hence, fluctuations in energy losses may be- 
come important due to this interaction. 



III. MONO-ENERGETIC TAU 

A. Propagated energy spectra 

Let us first consider the case of a mono-energetic beam 
of tau leptons propagating in rock. The simulated logio- 
distributions of the energy of an incident 3 • 10 9 GeV tau 
beam after crossing 1,5, and 10 km of standard rock are 
displayed in Fig. 0] (top). We show on the same plots the 
means of the energy distributions, which are the values 
used within the continuous energy losses approximation. 
The most probable values are always greater than the 
mean values of the distributions. After 1 km, a good 
part of the simulated events still carry a large fraction 
of the initial energy Eq, meaning that this fraction of 
particles did not undergo many interactions. However, 
the distribution is asymmetric and there is already a long 
tail of events undergoing hard losses. For larger paths in 
rock, fluctuations in the energy losses increase, resulting 
in a broadening of the distribution and a smoothening of 
the high-energy cutoff. 

The bottom panel of Fig. [4] shows the same quanti- 
ties for a higher initial energy of the incident r beam 
(Eq = 3 • 10 10 GeV). As expected from the discussion 
about the effective interaction lengths for y > y cu t, those 
distributions reflect the fact that the fluctuations are 
larger at higher energies due to the higher rate of pho- 
tonuclear interactions, which are precisely those leading 
to harder losses. 



B. Tau-lepton range 
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FIG. 4: Energy logio-distributions of taus of Eq =3- 10 9 GeV 
after propagation in standard rock for 3 different depths (1, 
5 and 10 km). Also shown are the mean values of the energy 
distributions. 



For a fixed initial energy of the tau, it is interesting to 
calculate the survival probability of the tau as a function 
of the traversed depth. In case of continuous energy loss, 
the survival probability is simply given by 

Psurv(E , x) = exp ( - / d " ) (8) 

V JO *dec{E u ) J 



with E u the solution of 

Jeo 1K E t) 

We also computed P surv for the case of stochastic losses. 
We plot in Fig. [5] the resulting curves, showing that in 
addition to broadening the distributions, the effect of the 
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FIG. 5: Survival probability in standard rock as a function of 
depth for different initial energies. Blue: continuous energy 
losses; Red: stochastic energy losses. 



fluctuations is globally to reduce the survival probabili- 
ties with increasing energy. Note also the presence of a 
small tail corresponding at higher ranges, which reflects 
the cases for which the number of interactions is lower. 
This decrease can be quantified through the range R(Eq) 
of the tau lepton which is simply : 

/>oo 

R(Eo) = dx P surv (E ,x) (10) 
Jo 

We show in Fig. [5] the corresponding range of the tau 
lepton in standard rock. Let us remind here that the 
weak interactions of the tau, which we have neglected 
thorough the whole paper, might change a little bit the 
picture at the highest energies. However, rather than giv- 
ing numbers to be taken at face values, we are interested 
here in analyzing the effect of the fluctuations, which ap- 
pear to slightly reduce the range of the tau as its energy 
increases. 



IV. POWER LAW INJECTION SPECTRA 

In this section, we consider the case of E~ x and E~ 2 
tau injection spectra. We restrict ourselves to these 
generic fluxes to study possible distorsions on the propa- 
gated fluxes due to the stochastic effects. All the calcula- 
tions we present here assume a continuous injection spec- 
trum following a power law between 10 8 GeV and 3 • 10 11 
GeV. The cutoff for the maximum energy is chosen to be 



FIG. 6: Range R(Eq) of the tau calculated with continuous 
(top curve) and with stochastic (bottom curve) energy losses. 
For comparison, the decay length of the tau (dashed line) is 
also shown. 



sharp to exhibit most clearly the different behaviors of 
the propagated fluxes. 

In Figs. [7] and [8j we show the results obtained respec- 
tively from E^ 1 and E~ 2 injection spectra. For a better 
lecture of the different structures, we choose to plot the 
fluxes corrected for the injection spectrum. This means 
that in a trivial situation with no losses and no decay 
of the tau, all results would be 1. Continuous lines are 
from the semi-analytical solution in the approximation of 
continuous energy losses, whereas squared dots are the 
results of the Monte-Carlo simulation including stochas- 
tic effects. Again, the sharp cutoff present in the semi- 
analytical solutions results from the choice of the maxi- 
mal energy together with the univocal relation between 
this maximal energy at x — and the propagated energy 
at a depth x. On the other hand, the stochastic treat- 
ment of the tau propagation has a smoothening effect on 
this cutoff: there are indeed fluctuations affecting a small 
fraction of particles which undergo less interactions and 
less hard losses. The energy range on which this broaden- 
ing occurs clearly increases with the depth traversed by 
the tau, and it is also more pronounced in the case of an 
E^ 1 flux. At the same time, we can see from both figures 
that this effect is the only important distorsion induced 
by the stochastic processes. This means that, as far we 
are dealing with continuous incident spectra behaving as 
power laws, there is a compensation between positive and 
negative fluctuations of the energy losses everywhere in 
the considered energy range, except near the high-energy 
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FIG. 7: Comparison of the semi-analytical solution (solid line) 
of the transport equation within the continuous energy losses 
approximation (Eqn[6]) with respect to a Monte-Carlo cal- 
culation (squared points) taking into account the stochastic 
effects of the radiative processes. An flux between 10 s 
GeV and 3 • 10 11 GeV is used at injection. Three distances of 
propagation in rock are shown: 1 km, 5 km and 10 km (from 
right to left). 



border. This compensation allows to use the continuous 
energy loss approximation as the correct mean value of 
the propagated spectrum. 



Energy [GeV] 

FIG. 8: Comparison of the semi-analytical solution (solid line) 
of the transport equation within the continuous energy losses 
approximation (Eqn[6j with respect to a Monte-Carlo cal- 
culation (squared points) taking into account the stochastic 
effects of the radiative processes. An E~ 2 flux between 10 s 
GeV and 3 • 10 11 GeV is used at injection. Three distances of 
propagation in rock are shown: 1 km, 5 km and 10 km (from 
right to left). 

other model leading to comparable or slightly greater rate 
would have to present a significantly harder differential 
cross section in term of y to challenge our conclusions. 



V. CONCLUSIONS 

In this paper, we have shown that the stochastic na- 
ture of the radiative processes undergone by tau leptons 
at ultra-high energy is indeed responsible for large fluc- 
tuations in the tau energy losses. At the same time, 
however, these fluctuations are not so large as to blur 
the picture with respect to the continuous energy loss 
approximation, as far as the calculation concerns power- 
law injection spectra in a given, continuous energy range. 

As already pointed out in Sec. Ill Ai significant theo- 
retical uncertainties exist in the calculations of the cross 
section for photonuclear interactions, which are the most 
relevant processes for tau energy losses at the ultra-high 
energies we are interested in. The present study was done 
on basis of one particular model [17| which is rather pes- 
simistic in the y range where the stochastic effects can 
lead to hard losses, in the sense that it leads to a fairly 
high rate of interactions. Any other model leading to a 
lower rate would not change the picture. Moreover, any 
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